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We study the effects of gauge-field fluctuations on the free energy of a homogeneous color super- 
conductor in the color-flavor-locked (CFL) phase. Gluonic fluctuations induce a strong first-order 
phase transition, in contrast to electronic superconductors where this transition is weakly first order. 
. . . The critical temperature for this transition is larger than the one corresponding to the diquark pair- 

, ing instability. The physical reason is that the gluonic Meissner masses suppress long-wavelength 

i fluctuations as compared to the normal conducting phase where gluons are massless, which stabilizes 

' the superconducting phase. In weak coupling, we analytically compute the temperatures associated 

with the limits of met ast ability of the normal and superconducting phases, as well as the latent 
O . heat associated with the first-order phase transition. We then extrapolate our results to intermedi- 

(]J ■ ate densities and numerically evaluate the temperature of the fluctuation-induced first-order phase 

I transition, as well as the discontinuity of the diquark condensate at the critical point. We find that 

, the London Umit of magnetic interactions is absent in color superconductivity. 

(N 



PACS numbers: 12.38.Mh, 24.85.+p 



I. INTRODUCTION 



(N 

■ It was suggested a long time ago that quark matter might exist within the central regions of superdense stars 

■ Since Quantum Chromodynamics (QCD) is an asymptotically free theory 2], it was argued that the extremely 
compressed matter found in neutron stars consists of quarks rather than of hadrons and that realistic calculations 

^ ■ in the framework of QCD become possible j^]. At high baryon densities and sufficiently low temperatures, however, 
Qh! a phase transition between normal and color-superconducting quark matter is expected 4, 5, 6, 7]. Therefore, 
color superconductivity (CSC) may be relevant to explain several important aspects of the highly compressed matter 
present in compact stars, e.g., the cooling rates Q, and the rotational properties of stars 0. Nevertheless, trustworthy 



. perturbative calculations can only be performed for ultra-high chemical potentials, /i ^ Aqcd- The weak-coupling 



expansion of the temperature for the diquark pairing instability reads 



al poter 



3 



TT 



where g is QCD running coupling constant at the chemical potential /i, and 7 is the Euler-Mascheroni constant. The 
energy gap at T = in the CFL phase is 

A(0) = 2-^7re-^Tc . (2) 

It follows from the Ginzburg-Landau (GL) theory of CSC in weak coupling that the phase transition is of second 
order at Tc and the GL parameter is ^12l] 

/ 727r3 Tc 
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with as = ^^/(47r). When ^ ^ then ^ and, therefore, the CFL color superconductor is of extreme type 1. 

It is weh known in the context of electronic superconductors that gauge-field fiuctuations change the second-order 
phase transition into a first-order transition However, the strength of this first-order transition is sensitive to 
the relationship among the three length scales that are involved: the coherence length near the transition, ^, the 
magnetic penetration depth near the transition. A, and the coherence length at T = 0, A superconductor 

with A ^ ^0 is said to be in the London limit. In this case, the coupling between the gauge field and the order 
parameter is approximately local. The opposite case, A <C ^o, corresponds to the Pippard limit and the coupling 
becomes highly nonlocal 14]. For a type-I electronic superconductor the Pippard limit is always realized at T = 0. 
As the temperature is raised towards the transition to normal quark matter, the penetration depth increases and so 
does the ratio A/^o- A crossover from the Pippard limit to the London limit would be expected if the transition is of 
second order. How does the first-order phase transition induced by gauge-field fiuctuations change this scenario? Are 
both limits still realized? In the case of known electronic superconductors of type I, the first-order phase transition 
is sufficiently weak to warrant a crossover between Pippard and London limits, which has been indeed observed 
experimentally for strong type-I materials like aluminum [Isl. However, the situation is completely different for a 
color superconductor. It was shown in a previous work 1^ that A <C is maintained at the phase transition for 
asymptotically high baryon density. As will be shown below, this feature remains valid when the results of Ref. 
are extrapolated to moderately high densities. 

The current work, which is an extension and continuation of the previous project is organized as follows. In 
the next section, we shall review the generalized GL free energy derived previously. The relevant thermodynamic 
quantities of the first-order color-superconducting transition will be calculated in weak coupling in Sec. IIIII and the 
extrapolation of the results to moderate coupling will be presented in Sec. II VI Concluding remarks will then be given 
in Sec.0 Moreover, technical details on the derivation of the generalized GL free energy, which were skipped in Ref. 
|l6( , will be sketched in the appendix. In contrast to Ref. |l6| , the zero-temperature coherence length will be defined 
as ^0 = l/(27rTc). Our units are h = c = ks = ^ and 4-vectors are denoted by capital letters, K = = (cj, k). In 
our formulas, Tr indicates the summation over all indices including momentum, /c, and energy, cj, while tr denotes 
the summation over all indices except momentum and energy. 



II. THE GENERALIZED GINZBURG-LANDAU FREE ENERGY 



The CJT effective potential 6] reads 

T[D,S] = ^ {Trlnp-^ + Tt{D-^V - 1) - Trln^"^ - TiiS'^S - 1) -2T2[D,S]} , (4) 

where Q denotes the 3-volume of the system, f> and S are the full gluon and quark propagators, and are 
the corresponding inverse tree- level propagators, and r2 is the sum of all two-particle irreducible vacuum diagrams. 
We work in the two- loop approximation, i.e., r2 contains only the diagrams shown in Fig. ^ The first diagram, 
containing quark propagators, leads to a term of order ^^/i^ in F, while the other two diagrams, containing only gluon 
propagators, lead to terms proportional to powers of T. Therefore, at small temperatures T ^ Tc ^ /iexp(— 1/^), one 
can drop the last two diagrams and restrict the consideration to the first. In explicit form, 

r2[A5] = -^Tr{Pn[5]}, (5) 

where 

H[5] = ]^Tr(VSVS) (6) 

is a functional of the full quark propagator S and the bare quark-gluon vertex F. Note that the trace in Eq. (|5J is 
over momenta, as well as Lorentz and adjoint color indices, while in Eq. it is over momenta, as well Dirac, fiavor, 
and fundamental color indices. The minus sign in Eq. Q takes account of the fermion loop and the factor 1/2 is due 
to the fact that this is a second-order correction to the CJT effective potential. The factor 1/2 in Eq. Q accounts 
for the doubling of the fermionic degrees of freedom in Nambu-Gor'kov space. 

The free energy is given by the CJT effective potential at its stationary points, determined by 
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FIG. 1: The two- loop approximation to Straight hnes denote quark propagators, wavy hnes denote gluon propagators. 

The first condition gives a Dyson- Schwinger equation for the gluon propagator, 

V-^ =D-^^U[S] . (8) 
One now sees that, at the stationary point, the second term in Eq. cancels the last term, 

r[P,5] = ^ [Trlnp-^ -Trln^-^ -Tr(5-^5- 1)] . (9) 

In terms of the gluon and quark propagators in the normal phase, Vn{K) and Sn{K)^ the propagators in the super- 
conducting phase are written as 

S{K) = Sn{K)^SS{K,A) , (10a) 
V-\K) = V-\K)^SU{K,A) , (10b) 

where SIl = 11 — 11^, i.e., already contains the hard-dense- loop (HDL) resummed gluon self-energy 11^. The 

gluon self-energy in the superconducting phase, 11, depends on the superconducting gap parameter. A, and so SU also 
depends on A. Similarly, the quark propagator in the normal phase Sn contains quark self-energy corrections, and 
6S depends on A. 

Inserting Eqs. (|TOj) into Eq. we obtain 

r = + Tcond + Tfluc + ^^Huc ' (11) 

where 

Tn = ^ [TVlnP-i - TVln5-i - TriS-'S^ - 1)] , (12a) 

Tcond = ^ [TV(P„OT)-TV(5-i(55)+Trln(l+ 5-1(55)] , (12b) 

Tfluc = ^ tr {In [1 + V„{K)5U{K, A)] - Vn{K)SU{K, A)} , (12c) 

k,uj=0 

Tfluc = ^ E [1 + ^niK)SIliK, A)] - V^iK)5UiK, A)} . (12d) 

The generalized GL free energy is the difference in the GJT effective potential between the superconducting phase 
and the normal phase, T — T^. It includes both the ordinary GL terms and the fluctuation terms. Note that we 
have added a term Tr(Pn^n) in T^nd and simultaneously subtracted it in rfluc,r|^^^. This term corresponds to the 
so-called exchange (free) energy jl7j and must be present in order to obtain the correct expression for Tcond !^ • 
Therefore, we have to subtract it in the fluctuation part of the free energy. Only with this subtraction, Ffiuc + ^^nuc 
represents the well-known plasmon ring resummation 17]. We note in passing that it is quite gratifying to see that 
the GJT formalism naturally contains all these different many-body contributions to the free energy. 

In Ref. j20|, the exchange energy was not subtracted from the plasmon ring contribution. This leads to an overall 
change of sign of the fluctuation energy. As shown below [see Eq. ([1^1) ]. the contribution of the fluctuation energy is 
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^ ln(l -^u) — which is always negative, while in Ref. it is ~ ln(l + u) which is positive (for u> 0). Therefore, 
the authors of Ref. 20] concluded that gauge-field fluctuations raise the free energy of the color-superconducting 
phase, and thus decrease the transition temperature to the normal phase. In our case, however, the gauge-fleld 
fluctuations decrease the free energy, i.e., stabilize the color-superconducting phase and therefore lead to a larger 
transition temperature. 

This is physically plausible if one remembers that gauge-fleld fluctuations are also present in the normal phase, 
namely in the flrst term in Eq. (|12a|) . Since transverse gluons are massless in the normal phase, 11^(0) = 0, long- 
wavelength fluctuations are enhanced over those in the color-superconducting phase where gluons are massive, SH ^ 0. 
Thus, the fluctuation energy in the normal phase is larger than in the superconducting phase. 

As shown in the appendix, the weak-coupling approximation gives rise to 



rcond 



21C(3) / M 



47r4 



f] A^(T), 



(13) 



where Tc is determined up to the accuracy of Eq. Q and t = {T — Tc)/Tc is the reduced temperature. The gap 
parameter of the fermionic quasiparticle excitations is A (8-fold) and 2 A (1-fold) 21]. One can check that the 
quadratic and quartic coefficients of Fcond for CSC are, respectively, 12 (= 8 x 1^ + 2^) and 24 (= 8 x 1^ + 2^) times 
larger than those for an electronic superconductor. The relevant fluctuation term is 



J w I 



In 



w?{T,k) 



w?{T,k) 



while rjjyj, is of higher order. The momentum-dependent Meissner mass reads 



(14) 



with the chromomagnetic penetration depth given by 

7C(3)^VA 
A2 247r4 T2 



(15) 



(16) 



and 



/(y) = 



7C(3) 



oo .1 



dx 



1 



■ X 



(5+|)[4(8+|)2+^2^21 



Carrying out the integration in Eq. (|14|) and combining the result with Eq. ([T^ . we flnd 

r2 



^iA2(T) + ^ 



A'^(T) + 327r(Te)'^F 



\\\T)) 



where the function F is deflned by 



F{z) 



dx X i In 1 



r/(^)} • 



The London limit corresponds to small arguments in Eqs. (|17|) and (|T^ . We have 



m = 1 



^C(5) 2 

140 C(3)^ 



0{y% 



and 



F{z) 



TT 3 



In the Pippard hmit the arguments of f{y) and F{z) become large and we end up with 

37r3 



28C(3)y 



^^41n, + 41n2+7^^^^_,^- 



(17) 



(18) 



(19) 



(20) 



(21) 



(22) 
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and 



Here we have retained the first corrections for both limits of the function f{y) in order to assess the deviation from 
each limit at the CSC phase transition. 

Before concluding this section, let us clarify once more several differences between our formulation and that of 
Ref. [23. First, since in their treatment the term — m^(T, in Eq. (|14|) . which arises from the subtraction of 
the exchange energy, is missing, their formal power-series expansion for the fluctuation energy in terms of A starts 
already at quadratic order. This then gives rise to a renormalized critical temperature T^. When including the term 
in question, there is no such renormalization of Tc- Nevertheless, the difference between Tc and is of order 0{g'^)^ 
and the two-loop approximation employed in the derivation of Eq. ([TH|) is not sufficiently accurate to provide all 
corrections of this order. Furthermore, the authors of Ref. |20| approximated the momentum-dependent Meissner 
mass by a constant and simply cut off the momentum integration in Eq. (|14|) . In that case, one is effectively in the 
London limit of Eqs. ([T7)) and ([TUI) . where the fluctuation energy is of the form of Eq. The shift in the transition 

temperature compared to that of Eq. (^) is then only of order 0{g^) 120|, and not of order O(^), as found here and 
in Ref. gl6i] . 

III. THE FIRST-ORDER CSC TRANSITION IN WEAK COUPLING 

A generic first-order phase transition can be described by three characteristic temperatures: the transition temper- 
ature T*, the maximum temperature of the (metastable) superheated superphase Tgh, and the minimum temperature 
of the (metastable) supercooled normal phase Tgc, respectively. These temperatures are related in the following way, 

< < Tsh , (24) 

and they can be obtained from the generalized GL free energy (|T^ . The lower margin of a supercooled normal phase 
corresponds to 



= 0, (25) 

A=0 



and, using Eq. H18|l . we have 

Tsc = Tc , (26) 
which relates Tgc with the onset temperature for diquark pairing. On the other hand, the transition occurs at 



OA 

for a value of A = Ac* 7^ 0. This implies that 



dT 

^ = , r = , (27) 



7C(3)A2 7C(3) 2w/eo' 



7C(3)Ag, 7C(3) /eg 



Eliminating t* in the equations above we have 



C§ \ 2167r^ /Te^' 



where J^{z) = -F'{z)/z + F{z)/z'^. Solving Eq. ^ for A^*, with the aid of Eq. (O, we obtain 
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The transition temperature is obtained substituting Eq. ((517)) into either one of Eqs. ((^ . which produces 

T; = (i + ^.)t.. (31) 
These were the results reported in Ref. . The penetration depth at the transition is 

which yields the ratio 



A?* 8647r4T2 



> 1 . (33) 



Thus, the Pippard limit is valid for the entire CSC phase at sufficiently large chemical potentials. 

We shall proceed to determine Tgh- The free energy F as a function of A has a local maximum between A = and 
the minimum Ac* at T = T* in the superconducting phase. As T increases, the local minimum remains unchanged 
until it coalesces with the local maximum, where 

dV d'^T 
for a value of A = Agh 7^ 0. It then follows that 



7C(3)A2, 7C(3) . 

+ + T8^^^^ ( — 1-0 



Moreover, Eq. ((S5|) together with Eq. yields 



= t:^^.9't'^' = Ul, . (36) 



126C(3)'^ " 2 

Subtracting the first equation in Eq. from Eq. IpT?^ and using ipHji . we find that 



^2 



U-t: = f^(l-ln2), (37) 



and as a result 

^sh 



l + ^(l-ln2) 



Note that T* is one order of g closer to Tgh than to T^c The ratio 

17287r4T2 



r; . (38) 



(39) 



implies that even the metastable CSC state is in the Pippard limit in weak coupling. Although the diagrammatics 
behind the generalized GL free energy function (|T^ determine Tc only up to subleading order, the leading-order 
differences among the three characteristic temperatures do not change if higher-order corrections to Tc are included. 

Another observable associated with the first-order phase transition is the latent heat L = TcAS, where AS is the 
change in entropy density at the transition. We have 



and as a result 

^'^ 2^(3)"^ " 7r2 



L = i^yn ^K^i'- (41) 
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FIG. 2: Comparison between the different temperatures involved in the discussion of the fluctuation-induced first-order phase 
transition. 



Now we calculate the strength of the first-order phase transition as was defined in Ref. IT^ . 

tuLM = 4^ , (42) 

where Ac^ is the jump in specific heat at the second-order phase transition, ignoring the fluctuations. If we ignore 
the third term in Eq. (|18|) we recover the ordinary GL theory from which we find Ac^ = 24/i^Tc/[7C(3)]. Thus, we 
have 

Note that Eq. (|43|) implies that the strength of the first-order phase transition weakens (logarithmically) with in- 
creasing chemical potential, which is in agreement with the fact that the second-order phase transition is recovered 
at asymptotically large densities. Note that for electronic superconductors, ^hlm/^c ^ 10~^ J^] which, for realistic 
values of ^ ~ 1 is much smaller than the right-hand side of Eq. ((I3j) . 



IV. NUMERICAL RESULTS 



Strictly speaking, the weak-coupling results in the previous section are only valid at ultra-high baryon densities 
such that ji ^ Aqcd- For quark matter that may exist inside a compact star ja is expected to be slightly higher 
than Aqcd, and then the weak-coupling expansion becomes problematic. Nevertheless, we shall assume that the 
generalized GL free energy remains numerically reliable down to realistic quark densities. Even though this is not the 
case, the qualitative statement for the absence of the London limit in CSC may still survive, due to the reason given 
at the end of this section. 

We solved Eq. (|29j) and the second equation in Eq. numerically in order to find Ac* and Agh as functions 
of the chemical potential. The transition temperature T* is obtained using Ac* in either one of Eqs. (|28|) and the 
temperature T^h is obtained from the first equation in Eq. We use the 3-loop formula for as = g'^/^Tr ^22] . 



8 



1.18 




o 



1.16 



1000 



le+06 



le+09 



le+12 1000 



le+06 



le+09 



le+12 



|Li (MeV) 

FIG. 3: (a) Comparison between the critical temperatures at high densities, (b) Discontinuity of the gap at the transition. 



47r 



/3oln(MVA|cD) 



2/3i ln[ln(^VA^cD) 
Pi ln(MVA|cD) 



/?4[ln(^VA^CD)]^ 



In 



In 



A2 



1/2 



(44) 



where f3o = 9, A = 51 - 19/(3A^/) = 32, 02 = 2857- 50337V//9 + 3257V|/27, for three colors and three flavors Nf = 3. 
Moreover, we have taken Aqcd = 364 MeV in our calculations, in order to obtain the correct value of ag at the scale 
of the Z-boson mass. 

Figure El shows the three temperature s TL c, T* and Tgh as functions of the chemical potential, along with the weak- 
coupling formula ((3T|) obtained in Ref. Note that T* is still closer to Tgh than to Tgc down to few hundreds of 
MeV. A comparison between the numerically evaluated critical temperature T* and Tc is shown in Fig. [HI (a) and the 
discontinuity of the gap at T*, relative to its value at T = 0, ((2j, is shown in Fig. [21(b) |23j. Both plots indicate that 



lim 



^ C 



= 1 



lim ——T 
M^oo A(0) 



= 0, 



(45) 



as expected, although the convergence is rather slow. 

Now we will address the question of whether or not the London limit is realized near T* for color superconductors 
in the range of chemical potentials explored here. From Fig. ^ we see that the ratio ^o/^c* ^ 1, meaning that only 
the Pippard limit of magnetic interactions is present in color superconductivity. Even for the minimum value of the 
ratio ^o/Ac*, which is around fi = 700 MeV, the Pippard expansion of m?{k^T) in Eq. works better than the 
London expansion, displayed in Eq. ((^ . This is also the case for the metastable CSC state up to Tgh, as is shown in 
Fig.ISl 

It is instructive to express the right-hand side of Eq. (|29|) in terms of the GL parameter and compare it with 
the corresponding equation for a metallic superconductor, whose generalized GL free energy was given in Ref. |T^ . 
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FIG. 4: (a) ^o/Xc* at small chemical potentials, (b) The same ratio at very large chemical potentials. 



We have 

for color superconductors and 



:f{ 



(a!:) " ' ^^^^ 



for electronic superconductivity. A large value on the right-hand side of Eq. or Eq. ((17)) points to the London 
limit at the first-order phase transition. Since Ofg <^ ag and Vi? ~ ag, the right-hand side of Eq. ((47|) is much larger 
than that of Eq. ((46|) under the same GL parameter. In other words, the London limit is more likely to be realized 
in metallic superconductors. 



V. CONCLUDING REMARKS 



In this paper we have systematically calculated the effects of gauge-field fluctuations on the free energy of a 
homogeneous CFL color superconductor in the Hartree-Fock approximation. We evaluated both analytically and 
numerically the temperature of the fluctuation-induced first-order phase transition, the latent heat, as well as the 
maximum temperature of a superheated superphase. It was also shown that the London limit for color-magnetic 
interactions in CFL color superconductors is absent. As the main reason we identified the weakness of electromagnetic 
interactions in comparison to strong interactions, ae <C ag. Thus, one can say that once the gauge-field fluctuations 
are taken into account, the local-coupling approximation between the order parameter and the gauge fields is not 
valid in the CFL phase. 

By using an inhomogeneous GL theory, lida and Baym 24] investigated the formation of vortices and supercurrents 
induced by external magnetic fields and rotation in pairing states near the critical temperature. Since they used a 
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FIG. 5: (a) ^o/Xsh at small chemical potentials, (b) The same ratio at very large chemical potentials. 



mean-field approximation, all gauge fields were regarded as averaged quantities and fluctuations around their mean 
values were not considered. In order to see how the inclusion of fluctuations would change their results one has to 
derive an effective action that depends only on the order parameter and the gauge fields. This action would display 
non-local interactions between the gau ge fi elds and the diquark condensate. Such an effective action could be obtained 
using the formalism developed in Ref. [25^ . 

It was shown in Ref. T^, by a one- loop renormalization-group calculation using the e expansion, that no stable 
infrared fixed point can exist for a theory involving local interactions between abelian gauge fields and order pa- 
rameters, unless the number of order parameter components. A/", is artificially extended to N > Nc = 365, which is 
far beyond the case of relevance for electronic superconductivity. This is then interpreted as signaling the presence 
of a first-order transition. Therefore, for electronic superconductors, gauge-field fluctuations are always expected to 
change the order of the phase transition to first order, irrespective of further details about the transition. For color 
superconductors the effective action containing only the order parameter and the gauge fields as well as the specific 
form of their interactions is not known, and the general result derived in Ref. may not be applicable. However, 
the results we obtained for the CFL phase seem to suggest that fluctuation-induced first-order phase transitions are 
indeed present in color superconductivity. Furthermore, due to the absence of the London limit, we expect that, once 
gauge- field fluctuations and first-order phase transitions are taken into account, local diquark-gluon interactions are 
never realized in color superconductors, regardless which phase is considered. This would constitute a striking new 
physical effect that would only come about in color superconductivity. In fact, the crossover from nonlocal to local 
interactions near the critical temperature in superconducting metals of strong type I has been recently observed . 
What we found in this paper rules out the possibility of observing such a crossover in color superconductors. 

Recently, a GL free energy that takes into account the effects of nonzero quark masses and charge neutrality has 
been derived within the mean- field approximation 26]. A study on the validity of local diquark-gluon interactions in 
this case and the effects of gauge-field fluctuations on the phase diagram obtained in Ref. |2^] is in progress and will 
be reported elsewhere p^ . 
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APPENDIX A 

In this appendix, we sketch some important steps for the derivation of the generalized GL free energy in the presence 
of gauge-field fiuct nations, which is shown in Eq. (|TKj) . in terms of Feynman diagrams. We have 

[S~\P)]Zl = {il^.P^.-liliP^W'5f,f,, (Ala) 
[5S-\P)]ll = i<P{Ph,p2{5}l6}l - 6%6fJ , (Alb) 

vUkU = 0, (Aid) 
vUkU - (Ale) 

where P = K = (cc;, ^), m|) = 3 ^^/i^/(27r^), p represents Pauli matrices with respect to Nambu-Gorkov indices, 

Ci and I' stand for the fundamental and adjoint color indices, respectively, fi are fundamental fiavor indices, and 
fi correspond to discrete Matsubara frequencies. The symbol in the gluon propagator means that we used the 
approximation for the total HDL gluon propagator that is relevant for the CSC energy scale. 

Diagrammatically, Vn is denoted by a wavy line, Sn is represented by a thick line, and the CSC correction to the 
inverse quark propagator (Alb) is associated with a two-point vertex bearing a cross. The corresponding diagrammatic 
expansions for 5S and JH are 



ss ^ + — ^ ^ 



sn 



In weak coupling we expand Sn as 

s- = 



12 



Expanding Fcond up to the fourth power of 0(P) we find 



= - 1/4 




+ 1/4 




+ 1/2 




1/2 




+ 3/8 




+ 3/2 



1/4 



where the weak-couphng approximation has been employed in order to retain the diagrams with at most one HDL 
gluon fine. The diagram bearing two crosses yields the expression ^pp,^{P)K[P\P')^{P'\ where the kernel K[P\P') 
is isomorphic to the kernel in the Dyson- Schwinger equation for the diquark scattering amplitude in the normal phase. 
Moreover, taking <^>(P) to be proportional to the pairing mode, i.e.. 



0(P) = A sin 



3^27 



■In 



(A2) 



where A is the energy gap, v = (3/2)^/^ ^^i^/(2567r'^/i), we have that 





2^ 



is proportional to T — Tc^ with Tc determined up to subleading order in g [see Eq. (^]. For the diagrams with four 
crosses the same mechanism yields 





at T = Tc, which reduces the number of quart ic terms in Fcond- Moreover, it will be shown at the end of this appendix 
that the following two diagrams 





are of higher order in weak coupling and can be dropped. For Fcond we end up with 



(A3) 



r = -1/4 

cond 




+ 1/4 




+ 1/2 




- 1/8 




which produces the terms in Eq. (^HJ- Now we treat the fiuctuation terms Ffiuc + F|^^^ , whose diagrammatic 
representation is 




- 1/3 /\ - 1/4 




(A4) 
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where Ffluc includes only the contribution from the static gluons and T'^^^ contains remaining contributions. Due to 
the Meissner effect, the shaded bubble does not vanish when the spatial momentum of the gluon line goes to zero at 
zero Matsubara energy. A resummation of all ring diagrams in Eq. (|A4|) is necessary for Ffluc and the result is the 
right-hand side of Eq. (|T^ . Regarding F^^^^, where the Matsubara energy of the gluon line is nonzero, dynamical 
screening prevents an infrared divergence for the integral over gluon momentum. In weak coupling T'^^^ is dominated 
by the first diagram in Eq. (|A4|) . which is again of higher order. Therefore, the contribution of F^^^ can be safely 
neglected. 

Now we present the argument supporting our assertion that the two diagrams in Eq. (|A3|) and the first diagram in 
F^^^ can be neglected in weak coupling. Let us denote the contribution of the first diagram in Eq. (|A3|) by CiA^. It 
consists of five free quark propagators with four- momentum {ly^pj and a self-energy insertion S(P) ~ g'^ u \n {jii / ~ 
gu. The main contribution to the ^-integration comes from a shell of thickness ~ \u\ around the Fermi surface and 
then we have 



^11 C 

which is of 0{g) in comparison to the quartic term in Eq. ([T^ . The contribution of the second diagram in Eq. ()A3|) . 
denoted by C2A^, can be estimated similarly. As is the case with the gap equation, the dominating contribution 
comes from the magnetic gluons with nonzero Matsubara energy. The integration for the quark propagators over the 
magnitude of their momenta p on each side of the gluon line can be approximately decoupled from the integration 
for the gluon propagator over the angle between p and p where the latter produces the forward logarithm. We then 
find 



c, ^ g'T'^' E ^ A 1^ f ^) - 3 ^, (A6) 

which is again of higher order. Now we consider the first diagram in F^^^ and denote its contribution as caA^. Since 

the typical momentum for the gluon line is k ^ m^^|cc;|^/^ co and ~ Tc, each bubble can be approximated by the 
static magnetic self-energy of gluons at the Pippard limit, i.e.. 




The sum over uj has a cutoff when uj ~ tud and then we end up with X^^-^qI^I ^ {l^l'^c) ^ 1/^- Consequently, 
we have cs ^ g fJ^'^/T^^ which is also negligible. 
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